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Based on an experimental study of two-body and three-body collisions in ultracold strontium 
samples, a novel optical-sympathetic cooling method in isotopic mixtures is demonstrated. Without 
evaporative cooling, a phase-space density of 6 x 10 -2 is obtained with a high spatial density that 
should allow to overcome the difficulties encountered so far to reach quantum degeneracy for Sr 
atoms. 
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The combination of laser cooling Q and evaporative 
cooling allowed to reach ultra-low temperatures and to 
observe Bose-Einstein condensation (BEC) and Fermi de- 
generacy for different atoms [2[. Laser cooling is indeed 
very effective to reach phase-space densities within few 
orders of magnitude of quantum degeneracy. The limits 
in cooling at high density are set by the optical depth of 
the sample, i.e. reabsorption of scattered light, and light- 
assisted atom- atom collisions. Forced evaporative cool- 
ing |3| represents the common way to circumvent these 
limits. However, this procedure does not work for all 
atoms. In particular, among the atoms cooled with opti- 
cal methods, none of the alkali-earth atoms reached quan- 
tum degeneracy so far, except ytterbium which has an 
alkali-earth-like electronic structure Hk A phase-space 
density of ~ 10 _1 was reported for Sr Q but BEC could 
not be reached. Sr is of interest for several reasons: In a 
natural strontium sample, different isotopes can be inves- 
tigated with different statistics; 88 Sr (82.6 %), 86 Sr (9.8 
%), and 84 Sr (0.6 %) are bosons with zero nuclear spin; 
87 Sr (7%) is a fermion with nuclear spin 1 = 9/2. A BEC 
of Sr atoms would allow the study of 0-spin condensates, 
fast optical cooling and continuous coherent matter- wave 
sources, quantum devices both for frequency standards 
and inertial sensors. Besides the studies towards quan- 
tum degeneracy, Sr has recently been the subject of ac- 
tive research in various fields spanning from laser cooling 
physics 0, 0, investigation of ultra-narrow transitions 
towards future optical clocks 13. @. ITol ITU . multiple scat- 
tering [13, and collisional physics [[I Hi [T^ [l7j. 

In this article, we present a novel scheme for the pro- 
duction of a high phase-space density sample of atomic 
strontium. Since collisional processes play a crucial role 
and relevant parameters are not well known for this atom, 
we studied 2-body and 3-body collisions both in pure 
samples, and isotopic mixtures. Our results show that 
usual evaporative cooling methods are difficult to ap- 
ply in a single species sample of 86 Sr or 88 Sr. The in- 
vestigation of an ultracold 86 Sr- 88 Sr mixture revealed a 
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large interspecies elastic cross-section; this allowed us to 
demonstrate an optical-sympathetic cooling method con- 
cept uallydifferent from what was done so far on neutral 
atoms pl| [Hi and somehow analogous to the one used 
for trapped ions [2p| , achieving a phase-space density of 
6xl0 -2 in conditions suitable to further cooling of the 
sample towards BEC. 

The experimental setup, based on the apparatus pre- 
viously described in 0, |2l| , allows to trap single-isotope 
samples or isotopic mixtures of 86 Sr and 88 Sr atoms in 
an optical dipole trap. The loading procedure starts with 
a magneto-optical trap on the 1 So- 1 Pi resonance transi- 
tion at 461 nm. During this phase, a decay channel from 
the excited state results in the accumulation of atoms 
in the magnetically trapped metastable 3 P2 state |2^ . 
The isotopic mixture is produced in the magnetic trap 
by stepping the frequency of the trapping laser to the ap- 
propriate values. After loading the magnetic trap, atoms 
are optically pumped back into the electronic ground 
state and a second cooling stage is applied based on a 
magneto-optical trap operating on the narrow 1 So- 3 Pi 
transition at 689 nm. At this point we typically obtain 
10 7 88 Sr atoms at 2 /iK, and 10 6 86 Sr atoms at 1 /iK. 
The isotopic mixture is then loaded into an optical dipole 
trap by superposing a focused 922 nm laser on the red 
magneto-optical trap, resulting in a 90 /iK deep poten- 
tial with radial and axial frequencies of 2 kHz and 26 
Hz, respectively. To optimize the transfer efficiency from 
the red magneto-optical trap into the optical dipole trap, 
the optical dipole beam waist is axially displaced from 
the magneto-optical trap by 500 jam (corresponding to 
2/3 of the dipole beam Rayleigh range and two times 
the red magneto-optical trap radii) in order to increase 
the magneto-optical trap-dipole trap spatial overlap [23| . 
The measured 6 s lifetime in the optical trap is consistent 
with the residual background gas pressure in our appara- 
tus. The off-resonant photon scattering induces a heating 
of 0.4 /iK/s. The polarization of the dipole beam is set to 
nearly fulfill the requirements for trapping at the " magic 
wavelength" that cancels the differential Stark shift for 
the 1 So- 3 Pi transition |24|. The atom number and tem- 
perature are measured independently on the two isotopes 
by absorption imaging with a resonant probe beam, and 
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the contribution of the non-resonant isotope is taken into 
account. The spatial density is inferred from the mea- 
sured temperature and the trapping potential. 

Two-body and three-body collisions are studied by 
loading single species or the isotopic mixture into the 
optical dipole trap. For the investigation of 2-body elas- 
tic collisions, the gas is initially put out of equilibrium by 
a sudden perturbation of the optical potential (pure sam- 
ples), or by an isotope-selective laser cooling pulse (mix- 
ture). We relate the measured thermalization time at the 
given density and temperature to the elastic cross-section 
g following and generalizing the approach described in 
[25^ . The resulting values are 088-88 = (3 =b 1) x 10 -17 
m 2 , (786-86 = (1.3 ± 0.5) x 10~ 14 m 2 , and a 8 6-88 = 
(4± 1) x 10 -16 m 2 . In the single isotope case, the results 
are consistent with the corresponding scatte ring lengths 
measured by photoassociation spectroscopy [T^J, [13 • No 
previous reference data is available for the two-isotope 
collisions. 

We studied the 3-body recombination rate K 3 by load- 
ing a single isotope (either 86 Sr or 88 Sr) into the dipole 
trap and observing the evolution of the atom number 
and temperature. With 88 Sr we found no experimental 
evidence for non-exponential decay and, given the initial 
atom density at the trap center of 3 x 10 13 cm -3 , we set 
an upper bound on if 88 at 10 -27 cm 6 /s. To determine 
iff 6 we followed the analysis described in 26]: the loss 
rate due to 3-body recombination is modeled by the rate 
equation 

d -^ = -K,J v n\r,t)d\ (1) 

or equivalently 

In N ® K C A*> [ ^Ml/3 r ( 9) 

ln m-- K3 Jo dt Jv^w dr ' (2) 

where n is the density, and N(t) is the atom number 
at time t. Figure [I] shows the data and the fit using 
equation[2|from which we obtain the rate constant K 86 = 
(1.0±0.5) x 10 -24 cm 6 / s. This recombination rate is more 
than three orders of magnitude larger than the typical 
values for ultracold alkali vapors J27J • 

The first consequence of these results is that evapora- 
tive cooling on pure samples of either 86 Sr or 88 Sr can- 
not be efficient. 86 Sr presents an extremely large elastic 
cross-section, a good point to establish a fast thermaliza- 
tion, but the 3-body recombination rate introduces a loss 
channel that is fatal with the typical geometries accessi- 
ble through optical dipole trapping. An optical trap with 
a much larger trapping volume would partially suppress 
this loss channel [2£|. 88 Sr instead turns out to be stable 
against 3-body decay, but the small elastic cross-section 
results in a long thermalization time compared to typi- 
cal trap lifetime. Indeed, in our experiment evaporative 
cooling of 88 Sr by reducing the dipole beam intensity al- 
lowed to increase the phase-space density to a maximum 
of 2 x 10 _1 , limited by the reduction of the thermalization 
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FIG. 1: Non-exponential decay of the number of 86 Sr, 
trapped in the optical dipole trap due to 3-body recombina- 
tion. Inset: the 3-body recombination rate constant is given 
by the slope of the natural log of the atom number with re- 
spect to jl(n 2 (t'))dt' where {n 2 (t')) = — ^ J v n 3 (f, t')dPr 
(see the text). The term t/r accounts for losses due to back- 
ground collisions. 



rate during the evaporation and by the 6 s background 
limited lifetime. 

These results led to a new all-optical sympathetic cool- 
ing scheme. The basic idea of the new cooling method is 
to use a mixture of isotopes trapped in the dipole trap: 
one isotope is laser cooled and by elastic collisions cools 
the other isotope which is the one of interest. The ad- 
vantage, compared to the schemes used so far for neutral 
atoms, is that the limitations of laser cooling of spatially 
dense and optically thick samples are drastically reduced. 
In this work, continuous laser cooling of 86 Sr leads to heat 
dissipation in a 88 Sr sample by sympathetic cooling; the 
small 86 Sr optical depth does not limit the achievable 
minimum temperature. Sympathetic cooling with neu- 
tral atoms normally requires a thermal bath with a large 
heat capacity with respect to that of the sample to be 
cooled. This is due to the fact that when the thermal 
bath is cooled by evaporative cooling, each lost atom 
carries an energy of the order of few times the temper- 
ature of the sample. In the case of optical-sympathetic 
cooling, each laser-cooled atom can subtract an energy 
of the order of the optical recoil in a time correspond- 
ing to a few lifetimes of the excited state (say 4, when 
operating at 1/4 I 8 at), without being lost. Considering 
the 1 So- 3 Pi transition, the power subtracted per atom 
is of the order of 6 mK/s; if the inter-isotope thermal- 
ization is fast enough, each 86 Sr atom can cool about 
10 3 88 Sr atoms down to /iK temperatures in one sec- 
ond. In the 86 Sr- 88 Sr isotope mixture the relatively large 
inter-species cross-section results in thermalization times 
typically of the order of few milliseconds. This thermal- 
ization is fast even on the time scale of laser cooling on 
the intercombination ^o- 3 ?! transition. On the other 
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FIG. 2: Dynamics of an optically trapped 88 Sr cloud sym- 
pathetically cooled with laser cooled 86 Sr. Filled circles •: 
88 Sr atom number. Open circles o: 86 Sr atom number. The 
number of 88 Sr atoms remains constant during the process, 
while 86 Sr decays exponentially with a 80 ms time constant. 
Under optimized conditions, the temperature (triangles V) 
decreases with a 150 ms time constant and the mixture is 
always at thermal equilibrium. 



hand, the 164 MHz 86 Sr- 88 Sr isotopic shift and the natu- 
ral linewidth of 7.6 kHz for the ^o^Pi transition, makes 
laser cooling on one isotope insensitive to the presence of 
the other. 

We implemented the optical-sympathetic cooling 
scheme by extending the temporal overlap between the 
optical dipole trap and the 86 Sr red magneto-optical trap 
after switching off the 88 Sr magneto-optical trap. During 
this phase the parameters that optimize the 88 Sr cooling 
are substantially the same as those typically employed 
in the second cooling stage (intensity 30-100 I sat with 
1 sat — 3/iW/cm 2 , -200 kHz frequency detuning, mag- 
netic field gradient 0.6 mT/cm). The horizontal dis- 
placement between the red magneto-optical trap and the 
dipole trap optimizes the spatial overlap and provides a 
continuous flux of 86 Sr atoms from the magneto-optical 
trap to the optical dipole trap replacing the atoms that 
are lost during the optical-sympathetic cooling due to 
light assisted collisions. Figure O reports the dynamics 
of optical-sympathetic cooling, starting from an initial 
temperature of 15-20 //K, limited by density dependent 
heating 0, [2<j . It can be observed that the cooling pro- 
cess does not induce losses on 88 Sr while the number of 
86 Sr atoms exponentially decays with a 80 ms lifetime. 
About 100 ms after switching off the 88 Sr red magneto- 
optical trap, we observe that the mixture attains thermal 
equilibrium. From this point on the equilibrium is main- 
tained, which is an essential condition for the implemen- 
tation of this cooling scheme. Under optimized condi- 
tions (overall optical intensity 100 I sat ) the temperature 
decays from the initial value with a 150 ms time constant. 
The minimum attainable temperature depends both on 
the intensity of the 86 Sr cooling beam, and the 88 Sr den- 
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FIG. 3: Asymptotic temperature of the 88 Sr cloud for differ- 
ent intensities of the 86 Sr laser cooling beam on the ^"So- 3 ?! 
transition (I sa t = 3 /iW/cm 2 ). The 88 Sr atom number is 
6 x 10 5 . 



sity. Figure 01 reports the asymptotic temperature as a 
function of the overall cooling light intensity incident on 
the sample, for a fixed amount of trapped 88 Sr. At high 
intensity, the 86 Sr lifetime increases due to the reduced 
spatial density and the final temperature is proportional 
to the intensity as expected from the usual Doppler the- 
ory. By lowering the intensity, the final temperature and 
the 86 Sr lifetime decrease to the point where 86 Sr is lost 
too fast with respect to the cooling dynamics of 88 Sr. 
The optimum values are obtained when the overall cool- 
ing light incident on the atoms has an intensity ~ 30 I sa t • 
For 6 x 10 5 88 Sr atoms trapped in the dipole trap, the 
final temperature is 6.7/iK at a peak density of 1.3 x 10 14 
cm 3 ; the corresponding phase-space density is 5 x 10 -2 . 
This value is only a factor of two lower with respect to 
what was previously obtained without forced evaporation 
|5| , but it exhibits more favorable conditions for starting 
evaporative cooling, considering the higher spatial den- 
sity (more than one order of magnitude higher) and the 
larger number of trapped atoms. By keeping the cooling 
parameters on 86 Sr fixed at the optimum value and by 
varying the number of trapped 88 Sr, we determined the 
dependence of the final temperature with respect to the 
88 Sr density. Figure^shows the dependence of 88 Sr tem- 
perature from the number of atoms in the trap. We deter- 
mine the density-dependent heating dT/dn ~2 /iK/(10 14 
cm -3 ), which is 20 times lower than the equivalent value 
for pure laser cooled 88 Sr [2^. This strong reduction 
is a direct consequence of the optical distinction of the 
two isotopes on the 1 So- 3 Pi transition. The limit on the 
86 Sr temperature of 4 jiK for zero 88 Sr density can be 
attributed to the laser cooling dynamics in the tightly 
confining potential of the dipole trap. Considering the 
power laws involved, reducing the optical dipole poten- 
tial in the final stage of the optical-sympathetic cooling 
should lead to further cooling of the sample and consid- 
erably increase the 88 Sr phase-space density. 

In conclusion, we investigated the collisional properties 
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FIG. 4: Temperature and phase-space density of the Sr 
cloud sympathetically cooled with laser cooled 86 Sr, as a func- 
tion of the 88 Sr atom number. 



measure a rather large inter-species elastic cross section. 
We exploit the fast 86 Sr- 88 Sr thermalization to demon- 
strate an optical-sympathetic cooling scheme in which 
heat dissipation for one isotope is obtained by laser cool- 
ing; the sympathetically cooled component is more than 
one order of magnitude more abundant than the laser 
cooled one, reducing substantially the limitations asso- 
ciated to laser cooling at high spatial density. We cool 
up to 7 x 10 5 88 Sr atoms at a peak spatial density of 
1.4 x 10 14 cm -3 and 6 x 10 -2 phase space density. The 
combination of this cooling scheme to a final evapora- 
tive cooling stage on pure 88 Sr should allow to increase 
considerably the phase-space density opening the pos- 
sibility to investigate quantum degenerate gases of Sr. 
An interesting possibility is also the development of new 
atomic sensors based on ultracold Sr atoms exploting the 
weak interatomic interaction and their small sensitivity 
to magnetic and electric fields. 



of ultracold Sr atoms of interest for the production of a 
Bose-Einstein condensate. On pure 86 Sr we find a very 
large 3-body recombination rate which sets stringent lim- 
its on the evaporation efficiency on this isotope. For 88 Sr 
we find that the small elastic cross-section results in a 
long thermalization time compared to typical trap life- 
time. This can explain why BEC of Sr was not achieved 
so far. On the other hand, for the 86 Sr- 88 Sr mixture we 
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